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Abstract 

A subset (~ 10%) of massive stars present strong, globally ordered (mostly dipolar) magnetic fields. The trapping 
and channeling of their stellar winds in closed magnetic loops leads to magnetically confined wind shocks (MOWS), with 
pre-shock flow speeds that are some fraction of the wind terminal speed. These shocks generate hot plasma, a source of 
X-rays. In the last decade, several developments took place, notably the determination of the hot plasma properties for 
a large sample of objects using XMM and Chandra, as well as fully self-consistent MHD modelling and the identification 
of shock retreat effects in weak winds. Despite a few exceptions, the combination of magnetic confinement, shock retreat 
and rotation effects seems to be able to account for X-ray emission in massive OB stars. Here we review these new 
observational and theoretical aspects of this X-ray emission and envisage some perspectives for the next generation of 
X-ray observatories. 
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1. Introduction 

Hot luminous, massive stars of spectral type O and B 
are prominent sources of X-rays which can originate from 
three distinct sources: shocks in their high-speed radia- 
tively driven stellar winds, wind-wind collisions in binary 
systems and magnetically confined wind shocks. 

In single, non-magnetic O stars, the intrinsic instabil¬ 
ity of wind driving by line-scattering leads to embedded 
wind shocks that are thought to be the source of their rel¬ 
atively soft (~0.5keV) X-ray spectrum, with a total X-ray 
luminosity that scales with stellar bolometric luminosity. 
Lx ~ 10“^ X Lboi (Chlebowski, 1989; Berghoefer et ah, 
1997; Naze et ah, 2011). In massive binary systems the 
collision of the two stellar winds at up to the wind ter¬ 
minal speeds can lead to even higher Lx, generally with a 
significantly harder (up to lOkeV) spectrum (for a review, 
see Rauw & Naze, this volume). 

Here we discuss a third source of X-rays from OB winds, 
namely those observed from the subset (~10%) of mas¬ 
sive stars with strong, globally ordered (often significantly 
dipolar) magnetic fields (Petit et al., 2013); in this case, 
the trapping and channeling of the stellar wind in closed 
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magnetic loops leads to magnetically confined wind shocks 
(MCWS) (Babel and Montmerle, 1997a, b), with pre-shock 
flow speeds that are some fraction of the wind terminal 
speed, resulting in intermediate energies for the shocks 
and associated X-rays (~2keV). A prototypical example 
is provided by the magnetic 0-type star 9^ Ori C, which 
shows moderately hard X-ray emission with a rotational 
phase variation that matches well the expectations of the 
MCWS paradigm (Gagne et ah, 2005a). 

Here, we discuss theoretical aspects of magnetic con¬ 
finement that determine the extent of the influence of the 
field over the wind. We, then, describe an effect called 
‘shock retreat’, which can moderate the strength of the X- 
rays, or even quench it altogether in extremely low mass 
loss rate stars, and the effects of rotation. In §3, we de¬ 
scribe current results in the X-ray observations of early- 
type magnetic stars, while §4 briefly discusses future out¬ 
looks of MCWS in X-ray astronomy. 

2. Theoretical Perspective 

To explain X-ray emission from the Ap/Bp star IQ 
Aur Babel and Montmerle (1997a) introduced the MCWS 
model. In their approach, they effectively prescribed a 
fixed magnetic field geometry to channel the wind outflow 
(see also Shore and Brown 1990). For large magnetic loops, 
wind material from opposite footpoints is accelerated to a 
substantial fraction of the wind terminal speed (i.e., >1000 
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km s“^) before the channeling toward the loop tops forces 
a collision with very strong shocks, thereby heating the gas 
to temperatures (10^ - 10® K) that are high enough to emit 
hard (few keV) X-rays. This star has a quite strong field 
(~4 kG) and a rather weak wind, with an estimated mass 
loss rate of about ~ yr“^, and thus indeed could 

be reasonably modeled within the framework of prescribed 
magnetic held geometry. However, the actual X-ray emis¬ 
sion from the star, which is predominantly soft, is further 
inhuenced by an effect called ‘shock retreat’ that we de¬ 
scribe below. Later, Babel and Montmerle (1997b) applied 
this model to explain the periodic variation of X-ray emis¬ 
sion of the 07 star 6^ Ori C, which has a lower magnetic 
held (^1100 G) and signihcantly stronger wind (mass-loss 
rate ~ IO^^Mq yr“^), raising now the possibility that the 
wind itself could inhuence the held geometry in a way that 
is not considered in the simple hxed-held approach. 


2.1. Magnetic Confinement 

In an interplay between magnetic held and stellar wind, 
the dominance of the held is determined by how strong it 
is relative to the wind. To understand the competition 
between these two, ud-Doula and Owocki (2002) dehned a 
characteristic parameter for the relative effectiveness of the 
magnetic helds in conhning and/or channeling the wind 
outhow. Specihcally, consider the ratio between the energy 
densities of held vs. how. 
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where the latitudinal variation of the surface held has the 
dipole form given by B^{6) = B^{cos^ 0-|-sin^ ^/4). In gen¬ 
eral, a magnetically channeled outhow will have a complex 
how geometry, but for convenience, the second equality in 
eqn. (1) simply characterizes the wind strength in terms of 
a spherically symmetric mass loss rate M = Airr^pv. The 
third equality likewise characterizes the radial variation of 
outhow velocity in terms of the phenomenological velocity 
law v{r) = Uoo(l — Rtf/rY, with v^o the wind terminal 
speed and /3 = 1; this equation furthermore models the 
magnetic held strength decline as a power-law in radius, 
B{r) = i3*(i?*/r)(”+^\ where, e.g., for a dipole n = 2. 

With the spatial variations of this energy ratio thus 
isolated within the right square bracket, we see that the left 
square bracket represents a dimensionless constant that 
characterizes the overall relative strength of held vs. wind. 
Evaluating this in the region of the magnetic equator (0 = 
90°), where the tendency toward a radial wind outhow is in 
most direct competition with the tendency for a horizontal 
orientation of the held, one can thus dehne an equatorial 
‘wind magnetic conhnement parameter’. 
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where M.g = M Mq/jt), B^oo = S „/(100 G), 
Ri 2 = i?*/(10^^ cm), and wg = Woo/(10® cm/s). In or¬ 
der to have any conhnement, p* > 1. As these stellar and 
wind parameters are scaled to typical values for an OB 
supergiant, e.g. C Pnp, the last equality in eqn. (2) im¬ 
mediately suggests that for such winds, signihcant mag¬ 
netic conhnement or channeling should require helds of 
order few hundred G. By contrast, in the case of the sun, 
the much weaker mass loss {Mq ~ 10“^^ Mq /yi) means 
that even a much weaker global held {Bo ~ 1 G) is suf- 
hcient to yield p^, ~ 40, implying a substantial magnetic 
conhnement of the solar coronal expansion. But in Bp 
stars the magnetic held strength can be of order kG with 
Mq ~ 10“^° Mq/jt leading ry* < 10®. Thus, the conhne¬ 
ment in Bp stars is very extreme. 

It should be emphasized that M used in the above for¬ 
malism is a value obtained for a spherically symmetric 
non-magnetic wind as the magnetic held may signihcantly 
inhuence the predicted circumstellar density and velocity 
structure. 


2.1.1. Alfven Radius 

The extent of the effectiveness of magnetic conhnement 
is set by the Alfven radius, Ra, where how and Alfven 
velocities are equal. This will also determine the extent 
of the largest loops and thus the highest shock velocities 
affecting the hardness of X-ray emission. This radius can 
be derived from eqn. ( 1 ) where the second square bracket 
factor shows the overall radial variation; n is the power-law 
exponent for radial decline of the assumed stellar held, e.g. 
n = 2 for a pure dipole, and with v{r) = '(;oo(l — R*/r)^ 
where f3 is the velocity-law index (typically /3 « 1). For a 
star with a non-zero held, we have 77 * > 0 , and so given 
the vanishing of the how speed at the atmospheric wind 
base, this energy ratio always starts as a large number 
near the stellar surface, p{r —>• i?*) —>• 00 . But from there 
outward it declines quite steeply, asymptotically as r~'^ 
for a dipole, crossing unity at the Alfven radius dehned 
implicitly by p{Ra) = 1- 

For a canonical /? = 1 wind velocity law, explicit so¬ 
lution for i?A along the magnetic equator requires hnding 
the appropriate root of 



(3) 


which for integer 2n is just a simple polynomial, specih¬ 
cally a quadratic, cubic, or quartic for n = 1, 1.5, or 2. 
Even for non-integer values of 2n, the relevant solutions 
can be approximated (via numerical htting) to within a 
few percent by the simple general expression, 

+ 1/4)1/(2») _ (l/4)l/(2n) ^ (4) 


For weak conhnement, ry* <C 1, we hnd Ra —> i?*, while for 
strong conhnement, 7y* ^ 1, we obtain Ra —t 
In particular, for the standard dipole case with n = 2, we 
expect the strong-conhnement scaling Ra/R* « pV^. 
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Clearly i?A represents the radius at which the wind 
speed V exceeds the local Alfven speed Va- It also charac¬ 
terizes the maximum radius where the magnetic field still 
dominates over the wind. For Ap/Bp stars where stellar 
fields are of order kG, r?* » 1, e.g. for a Ori E it is about 
10^, implying an Alfven radius ~ 60i?*. Thus, in Bp (and 
Ap) stars wind is trapped to large radii creating extensive 
magnetospheres. This also implies that X-rays from Bp 
stars should be intrinsically hard. But as we show in §2.2, 
shock retreat effects may soften it significantly. 

2.1.2. MHD Simulations 

The initial magnetohydrodynamic (MHD) simulations 
by ud-Doula and Owocki (2002) assumed, for simplicity, 
that radiative heating and cooling would keep the wind 
outflow nearly isothermal at roughly the stellar effective 
temperature. The simulations studied the dynamical com¬ 
petition between field and wind by evolving MHD simu¬ 
lations from an initial condition when a dipole magnetic 
field is suddenly introduced into a previously relaxed, one¬ 
dimensional spherically symmetric wind. 

Immediately after the introduction of the field, the dy¬ 
namic interplay between the wind and the field leads to two 
distinct regions. Along the polar region, the wind freely 
streams radially outward, stretching the field lines into a 
radial configuration, as can be inferred from the left panel 
of illustrative Fig. 1. If the field is strong enough, around 
the magnetic equator a region of closed magnetic loops is 
formed wherein the flow from opposite hemispheres col¬ 
lides to make strong shocks, quite similar to what was 
predicted in the semi-analytic, fixed-field models of Babel 
and Montmerle (1997a). The shocked material forms a 
dense disk-like structure which is opaque to line-driving. 
But, its support against gravity by the magnetic tension 
along the convex field lines is inherently unstable, leading 
to a complex pattern of fall back along the loop lines down 
to the star, again as suggested by the left panel of Fig. 1 . 

Note that even for weak field models with moderately 
small confinement, r;* < 1/10, the field still has a notice¬ 
able global influence on the wind, enhancing the density 
and decreasing the flow speed near the magnetic equator. 
However, shock speeds are probably not sufficiently large 
enough to produce any X-rays in that case. 

But to model the actual X-ray emission from shocks 
that form from the magnetic channeling and confinement, 
subsequent efforts (ud-Doula, 2003; Gagne et ah, 2005a) 
have relaxed the assumption of isothermal equation of state 
in earlier studies to include a detailed energy equation 
that follows the radiative cooling of shock-heated mate¬ 
rial. The magnetically channeled wind shock model pro¬ 
vides excellent agreement with the diagnostics from the 
phase-resolved Chandra spectroscopy of 0^ Ori C (see §3 
for more details). 

2.1.3. 3D MHD 

To get a more complete picture of the wind-field inter¬ 
action, one needs ultimately full 3D MHD simulations. Us- 
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Figure 2: Equatorial mass distribution drug/dr (in units Mq//?*) for 
the azimuthally averaged 3D model (left) and for a corresponding 2D 
model (right), plotted versus radius and time, with the dashed hor¬ 
izontal line showing the Alfven radius Rx ~ 2.23i2*. Following a 
similar adjustment to the initial condition, the long-term evolution 
of the 2D models is characterized by a complex pattern of repeated 
strings of infall, whereas the 3D azimuthally averaged model settles 
into a relatively smooth asymptotic state characterized by an en¬ 
hanced mass near and below the Alfven radius. Figure reproduced 
from ud-Doula et al. (2013). 

ing the specific parameters chosen to represent the proto¬ 
typical slowly rotating magnetic O star 0^ Ori G, for which 
centrifugal and other dynamical effects of rotation are neg¬ 
ligible, ud-Doula et al. (2013) have computed the first fully 
3D MHD model of its wind. The computed global struc¬ 
ture in latitude and radius resembles that found in previ¬ 
ous 2D simulations, with unimpeded outflow along open 
field lines near the magnetic poles, and a complex equa¬ 
torial belt of inner wind trapping by closed loops near 
the stellar surface, giving way to outflow above the Alfven 
radius. In contrast to this previous 2D work, the 3D sim¬ 
ulation also shows how this complex structure fragments 
in azimuth, forming distinct clumps of closed loop infall 
within the Alfven radius, transitioning in the outer wind 
to radial spokes of enhanced density with characteristic az¬ 
imuthal separation of 15° —20°. Applying these results in a 
3D code for line radiative transfer, they show that emission 
from the associated 3D ‘dynamical magnetosphere’ (DM) 
matches well the observed Ha emission seen from 0^ Ori G, 
fitting both its dynamic spectrum over rotational phase 
and the observed level of cycle-to-cycle stochastic varia¬ 
tion. Comparison with previously developed 2D models 
for the Balmer emission from a dynamical magnetosphere 
generally confirms that time averaging over 2D snapshots 
can be a good proxy for the spatial averaging over 3D az¬ 
imuthal wind structure, as illustrated in Fig. 2. Neverthe¬ 
less, fully 3D simulations will still be needed to model the 
emission from magnetospheres with non-dipole field com¬ 
ponents such as T Sco. A similar study examining X-rays 
predicted by this 3D MHD model is currently under way 
(Cohen et al., in preparation). 
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Figure 1: Colour plots of log density (left) and log temperature (middle) for arbitrary snapshot of structure in the standard model with 
r;* = 100. The right panel plots the proxy X-ray emission XEMt^ (weighted by the radius r), on a linear scale for a threshold X-ray 
temperature Tx = 1.5 MK. Figure reproduced from ud-Doula et al. (2014). 
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Figure 3: Schematic illustration of the “shock retreat” from ineffi¬ 
cient cooling associated with a lower mass loss rate M, showing a 
hemispheric, planar slice of a stellar dipole magnetic field. Wind 
outflow driven from opposite foot-points of closed magnetic loops is 
channeled into a collision near the loop top, forming magnetically 
confined wind shocks (MOWS). For the high M case in the upper 
panel, the efficient cooling keeps the shock-heated gas within a nar¬ 
row cooling layer, allowing the pre-shock wind to accelerate to a 
high speed and so produce strong shocks with strong, relatively hard 
X-ray emission. For the low M case in the lower panel, the ineffi¬ 
cient cooling forces a shock retreat down to lower radii with slower 
pre-shock wind, leading to weaker shocks with weaker, softer X-ray 
emission. Figure reproduced from ud-Doula et al. (2014). 


2.2. Shock Retreat 

Although magnetic confinement is the dominant effect 
in production of X-rays in magnetic massive stars, there 
are other secondary effects that can influence both the 
total luminosity and hardness of X-rays. In particular, 
an effect called ‘shock retreat’ directly related to cooling 
efficiency can significantly moderate X-rays in stars with 
low mass loss rates such as Bp stars. 

As illustrated schematically in Figure 3 (ud-Doula et al., 
2014, see also Figure 13 of Babel and Montmerle 1997a), 
the much lower mass loss rates of such B-stars also implies 
much less efficient cooling of the post-shock flow. When 
the associated cooling length, becomes comparable to the 
Alfven radius, the shock location is effectively forced to 
“retreat” back down the loop, to a lower radius where the 
lower wind speed yields a weaker shock, implying then 
a much softer X-ray spectrum. It was shown that such 
shock retreat can be effectively parametrized by a cooling 
parameter, Xoo oc V^/M. A high value of this cooling 
parameter, Xoo ■, implies an inefficient cooling, and, thus, a 
large shock retreat effect (ud-Doula et al., 2014). 

In extreme low mass loss cases, value of Xoo ^ 1 which 
may lead the shock retreat to be as large as the loop size. 
This can quench the wind within the closed magnetosphere 
suppressing the X-ray production altogether. Here, it may 
be interesting to note that confined winds were also sug¬ 
gested to play a role in the X-ray emission of mid-B to 
mid-A stars, when detected (Czesla and Schmitt, 2007; 
Stelzer et al., 2011; Robrade and Schmitt, 2011). Indeed, 
after examining the ratio of X-ray to radio luminosities, 
Robrade and Schmitt (2011) excluded low-mass compan¬ 
ions as the sources of X-rays. Considering a sample of 
Ap/Bp stars, these authors (see also Robrade 2014) fur¬ 
ther found that only stars brighter than ^ 2 OOL 0 emit 
X-rays, implying the presence of a threshold in physical 
phenomena, whose origin could be linked to shock retreat 
effects. 

Now, to estimate the total X-ray luminosity in typical 
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OB stars, Lx, ud-Doula et al. (2014) developed a semi- 
analytic analysis to derive a generalized “X-ray Analytic 
Dynamical Magnetosphere” (XADM) scaling law for X- 
rays emitted from such MOWS. To describe the paradigm 
briefly without much technical details, XADM assumes lo¬ 
cal latitudinal variation of radial mass flux at the stellar 
surface as prescribed by Owocki and ud-Doula (2004) . The 
size of the magnetosphere can be readily estimated from 
the magnetic confinement parameter, 77 *. It then computes 
shock velocity along a loop. Energy emitted as X-rays 
above a threshold can be estimated from kinetic energy 
dissipated in the shocks. Integration over the closed mag¬ 
netosphere then yields the total X-ray luminosity above a 
threshold assuming 100 % efficiency. 

The resultant XADM scaling follows a very similar 
trend to the full magnetohydrodynamic (MHD) simula¬ 
tion results, but is about a factor 5 higher. Compared to 
the idealized steady-state emission of the analytic XADM 
model, the numerical MHD simulations show extensive 
time variability with repeated intervals of infall of cooled, 
trapped material, and it appears that this lowers the over¬ 
all efficiency of X-ray emission to about 20 per cent of the 
idealized XADM prediction. 

2.3. The Effects of Rotation 

A further complication arises from rotation (not ac¬ 
counted for by XADM), as short-period modulations and 
substantial broadening in photospheric spectral lines (Conti 
and Ebbets, 1977; Fukuda, 1982), corresponding to pro¬ 
jected rotation speeds of hundreds of km/s, have been de¬ 
tected in some massive stars. 

For a non-magnetic rotating star, conservation of an¬ 
gular momentum in a wind outflow causes the azimuthal 
speed near the equator to decline outward as ^ 1 /r, 
meaning that rotation effects tend to be of diminishing 
importance in the outer wind. However, at the stellar sur¬ 
face, reduced effective gravity leads to enhanced mass loss 
rate (ud-Doula et al., 2008, 2009) and slower wind. 

By contrast, in a rotating star with a sufficiently strong 
magnetic field, magnetic torques on the wind can spin it 
up; for some region near the star, i.e., up to about the 
maximum loop closure radius i?c (closely related to Alfven 
radius, Ra), they can even maintain a nearly rigid-body 
rotation. As such, even for a star with surface rotation be¬ 
low the orbital speed, maintaining rigid rotation will even¬ 
tually lead to a balance between the outward centrifugal 
force from rotation and the inward force of gravity. 

Unsupported material at radii r < i?K, Kepler radius, 
will tend to fall back toward the star, but any material 
maintained in rigid-rotation to radii r > Rk will have a 
centrifugal force that exceeds gravity, and so will tend to 
be propelled further outward, unless held down by strong 
enough magnetic field as is the case for Bp and some other 
magnetic B stars. 

This complex interaction between the field and the 
wind naturally leads to two distinct groups of magnetic 


stars: “dynamical magnetospheres” (DM) and “centrifu¬ 
gal magnetospheres” (CM) (Sundqvist et al., 2012; Petit 
et al., 2013). For the former case, Ra < Rk and rota¬ 
tion is dynamically unimportant with all the supported 
material within the magnetosphere falling back toward the 
star. Notably, all Of?p stars fall within this category (their 
slow rotation probably comes from magnetic braking, see 
more details below). X-ray emission from DM is governed 
mainly by the magnetic confinement and shock retreat ef¬ 
fects. 

Whereas, for the latter case, Ra > Rk leading to co¬ 
rotating magnetospheres, as is the case for the rapidly ro¬ 
tating strong magnetic star a Ori E. The effect of rapid 
rotation on X-ray emission is currently a very much work 
in progress. However, preliminary studies show that the 
enhanced mass loss due to rotation leads to a more effi¬ 
cient cooling, which in turn can both increase and harden 
X-rays (Bard and Townsend, 2015). Here, there is a three 
way competition between rotation, magnetic confinement 
and shock retreat. Since the cooling parameter Xoo oc 
^shock/^’ the shock retreat effects might be moderated 
by the effects of rotation which increases M in general. 
But which effect ultimately dominates may depend on the 
details of stellar parameters of specific stars and rapidly 
rotating B stars can have a variety of X-ray luminosities. 
This is discussed further in §3. 

Another natural consequence of interaction of magnetic 
field with wind and rotation is the increase in angular mo¬ 
mentum loss from the stellar surface. It turns out that the 
angular momentum is carried away by not only the gas, 
but also by the field itself in the form of Maxwellian stress 
tensor. In a classical analysis, Weber and Davis (1967) 
showed that the total angular momentum loss can be ex¬ 
pressed as J = (2/3)MDi?^, where M is the mass loss 
rate, H is the stellar angular velocity, and Ra is a charac¬ 
teristic Alfven radius. ud-Doula et al. (2009) have shown 
that the typical timescale for spindown for massive mag¬ 
netic O stars is about 1 Myr or so. This has a direct effect 
on the evolution of these stars and their X-ray luminosi¬ 
ties. Currently, all known strongly magnetic 0-stars are 
slow rotators, with the exception of Plaskett’s star where 
interactions within the binary are speculated to play an 
additional role in spinup of the primary (Grunhut et al., 
2013; Palate and Rauw, 2014). 

3. Observational Perspective 

From the theoretical perspective, magnetic confinement, 
shock retreat and rotation effects are the main ingredients 
for production and explanation of X-rays in MOWS. Of 
course, other factors like the tilt angle of the field to the 
rotation axis or more generally the field geometry can also 
have substantial influence, as could also the presence of 
a companion, of stellar pulsations or of additional (az¬ 
imuthal) structures in the wind. Now, let us turn our 
attention to what the observations actually show. We 
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shall see that, despite a few exceptions, indeed these three 
parameters can mostly account for X-ray emissions from 
magnetic massive stars. 

The observed scatter in X-ray luminosity of AB stars 
may, as in OB stars, be related to different confinement, 
rotation, and cooling parameter values, but better knowl¬ 
edge of these parameters are required before a comprehen¬ 
sive comparison with models can be done just like for OB 
stars, so we will not describe these objects in detail here. 

3.1. X-ray luminosity 

The X-ray luminosity of massive OB stars is well known 


consistent with MHD models. The observed X-ray lumi¬ 
nosities for most of the known magnetic stars fall within 
this range. 

However, there are several clear exceptions (Fig. 4), 
five of which show observed luminosities larger than pre¬ 
dicted. A first outlier is r Sco (#11) which has a com¬ 
plex non-dipolar field, something the current version of 
XADM does not consider as it is limited to pure dipolar 
cases. Three cases, cr Ori E (#31), HD 182180 (#45) and 
HD 142184 (#47), are rapid rotators (Petit et ah, 2013). 
However, XADM does not take into account rapid rota¬ 
tion, but Bard and Townsend (2015) have shown that ro- 
to scale linearly with the bolometric luminosity, log[Lx/ABOL]"^^^^°^ both enhances and hardens the X-ray emission. It 


—7. However, there are some exceptions that display over¬ 
luminosities, like the well-known cases of 6^ Ori C or r Sco 
(Schulz et ah, 2000; Gagne et ah, 2005a; Mewe et ah, 2003; 
Cohen et ah, 2003). The discovery of magnetic fields in 
these objects is relatively recent, but they were found to 
play a key role in explaining their X-ray spectra. This 
overluminosity, of about one dex compared to the major¬ 
ity (“normal”) of OB stars, was further found in many 
other magnetic massive stars (u Ori E - Skinner et al. 2008 
or the Of?p stars - Naze et al. 2004, 2007, 2008, 2015; Petit 
et al. 2015). However, some other magnetic massive stars 
appear much fainter (Favata et al., 2009; Oskinova et al., 
2011; Ignace et al., 2013; Naze et al., 2014). 

Several attempts have been made recently to under¬ 
stand such diversity. A first, limited survey of 11 objects 
by Oskinova et al. (2011) found no link between the level 
of X-ray luminosity and the magnetic field strength, ro¬ 
tational period, bolometric luminosity, or pulsational pe¬ 
riod (when existing). The general trend amongst stars 
(more/less luminous) seemed however qualitatively accounted 
for by the formula of Babel and Montmerle (1997b), Lx oc 
MFooHO-A 

A larger survey, covering nearly two thirds of the cata¬ 
log of known magnetic OB stars of Petit et al. (2013) was 
undertaken by Naze et al. (2014). The 40 targets used in 
the survey were scattered all over the so called “magnetic 
confinement-rotation diagram” (see Petit et al., 2013), lim¬ 
iting the biases. The situation then appeared in more de¬ 
tail: some of the sources (0-stars and a few B-stars) have 
log[Lx/LBOL]~ —6.2 or equivalently Lx oc ® whereas 
other B-stars display Lx oc Lg^j or Lx oc (Fig. 4). 

These observations are clearly at odds with the earlier 
Babel & Montmerle prediction of Lx oc M. Moreover, 
the level of emission predicted by this formalism is 1.8 dex 
too high. However, the more recent XADM formalism, 
confirmed by MHD simulations, has shown much more 
promise (ud-Doula et al., 2014). The right panel of figure 4 
shows a comparison between the observed X-ray luminos¬ 
ity of magnetic stars and the predicted values using the 
XADM model of ud-Doula et al. (2014). The dotted line 
illustrates the ideal model with 100% efficiency whereas 
the solid line indicates a scaling by 10%; the grey shaded 
area corresponds to scalings by 5-20%, a range in efficiency 


may be noted in this context that those three objects are 
also particularly hard (Naze et al., 2014), but it should 
also be mentioned that the difference between rotating and 
non-rotating cases in Bard & Townsend models is less than 
one dex, while the observed difference for those three ob¬ 
jects is larger (about 2 dex). In addition, HD64740, a 
star with similar properties (rotation rate, confinement, 
cooling) as a Ori E does not display a bright and hard 
emission, and its Lx appears in line with XADM predic¬ 
tions. Finally, a last exception, HD200775 (#42) remains 
unexplained. 

There is also one case where the observed luminosity 
is smaller than the predicted one. It corresponds to NGC 
1624-2 (#5), which hosts a very strong (20kG) magnetic 
field. It turns out that the huge magnetosphere of this ex¬ 
treme magnetic star heavily attenuates its X-ray emission, 
making it appear less luminous and much harder than it in¬ 
trinsically is (Petit et al., 2015). In fact, the estimated in¬ 
trinsic luminosity is quite consistent with prediction from 
the XADM formalism. 

3.2. Absorption and hardness of spectrum 

Traditionally, a basic “magnetic-thus-hard” paradigm 
was assumed. However, one could a priori consider ob¬ 
served hardness of MOWS with care. Indeed, for NGG 
1624-2, the apparent hardness of the X-ray spectrum is 
actually due to the large absorption by the dense circum- 
stellar material (A# ~ 1 — 4 x 10^^cm“^), the intrinsic 
hardness being drastically different (Petit et al., 2015). 
MHD simulations of such an extreme magnetosphere are 
presently out-of-reach, but the general properties of this 
object qualitatively agree with those expected from con¬ 
fined winds in a extreme magnetosphere as predicted by 
XADM paradigm (Petit et al., 2015). 

However, NGG 1624-2 clearly is an extreme case. In the 
global survey of magnetic OB stars by Naze et al. (2014), 
the amount of local absorption in spectral fits was com¬ 
pared to that found from general surveys of OB stars and 
they were found to be similar. This suggests no signifi¬ 
cantly increased absorption due to confined winds in most 
magnetic OB stars, which is corroborated by the monitor¬ 
ing of 9^ OriG, where only a small increase (from 4.5 to 
5 X 10^^ cm“^, not even formally significant considering the 
errors) was found (Gagne et al., 2005a). This is consistent 
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Figure 4: Observed X-ray luminosity (corrected for ISM absorption) as a function of mass-loss rate (left) and of the predicted values using 
the XADM model of ud-Doula et al. (2014) (right). Left: Filled blue dots correspond to O-stars, black empty triangles to B-stars, magenta 
crosses and downward-pointing arrows to faint detections and upper limits on the X-ray luminosity, respectively. Blue and magenta lines 
represent Lx oc and Lx oc relation. Right: The dotted line illustrates the ideal model with 100% efficiency whereas the solid line 

indicates a scaling by 10%; the grey shaded area corresponds to scalings by 5-20% (a range in efficiency consistent with MHD models). The 
symbol shapes are the same as in the left panel, but this time symbols are colour-coded according to their measured hardness ratios. Figures 
reproduced from Naze et al. (2014), numbers correspond to star IDs defined in Petit et al. (2013). 



with the small density columns predicted by MHD models 
(see e.g. Fig. 10 of Petit et al., 2015). The observed hard¬ 
ness, once corrected for ISM absorption, is thus generally 
a good indicator of the intrinsic hardness. 

Observationally, the spectrum of 0^ Ori C (#3 in Fig. 4) 
is dominated by plasma at ~2.5keV (Schulz et al., 2000; 
Gagne et al., 2005a), and hot plasma is certainly present 
in other stars, e.g. rSco (Mewe et al., 2003; Cohen et al., 
2003) or Of?p stars (Naze et al., 2004, 2007, 2008). But 
the hard component does not always dominate, as in Of?p 
stars (which usually appear soft), or may not even be de¬ 
tectable (e.g. /3Cep, Favata et al., 2009). Even cases 
of “twins” wherein a pair of stars have twin-like similar 
stellar parameters, revealed some surprises, for instance 
with one r Sco analog being quite similar to its prototype 
(HD 63425) while the other one (HD 66665) was softer (Ig- 
nace et al., 2013). 

Surveys by Oskinova et al. (11 stars, 2011) and Naze 
et al. (40 stars, 2014) have noted such diversity of the 
situation. Comparison with general OB stars surveys fur¬ 
ther reveals that magnetic O-stars appear at least some¬ 
what harder than “normal” O-stars while magnetic B-stars 
seemed rather softer than usual B-stars (Naze et al., 2014). 
The latter conclusion may however be biased by the dif¬ 
ferent population sampled, the observed magnetic B-stars 
sampling lower X-ray luminosities than in surveys. Nev¬ 
ertheless, no correlation could be found between hardness 
and magnetic or stellar parameters (Naze et al., 2014). In 
particular, the prediction of harder X-ray emission in case 


of stronger magnetic confinement, larger mass-loss rate, or 
higher rotation rate was not detected. 

3.3. X-ray Variability 

By design, theoretical predictions such as XADM are 
time independent, although fully self-consistent MHD mod¬ 
els do suggest certain level of time variability in X-rays 
and other bands. Observationally, several types of vari¬ 
abilities associated with magnetically confined winds have 
been detected. In general, these variabilities do not change 
global properties of X-ray emission discussed above, but 
they do show how complex the interplay between the field 
and wind can be. 

3.3.1. Flares 

The first type of possible variations is short-term X-ray 
flaring. MHD models of dynamical magnetospheres show 
a complex temporal behaviour, with stochastic episodes of 
heating and infall back on the star - but they are limited 
to small regions, hence are smoothed out in 3D magneto¬ 
spheres (ud-Doula et al., 2013). In centrifugal magneto¬ 
spheres, however, the trapped material accumulates and, 
to keep the balance, some ejection needs to take place, 
through “centrifugal” breakouts involving magnetic recon¬ 
nection on the outer edges of magnetospheres (Townsend 
and Owocki, 2005; ud-Doula et al., 2006). This should gen¬ 
erate flares and such events have been recorded in cr Ori E 
by ROSAT, XMM, and Chandra (Groote and Schmitt, 
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2004; Sanz-Forcada et al., 2004; Skinner et al., 2008; Ca¬ 
ballero et al., 2009), but their interpretation was much 
debated: is it intrinsic to tr OriE or due to an unseen low- 
mass companion? Several arguments can be put forward 
for the former scenario. For example, the decay part of 
the flare lightcurve suggests a size of ~ li?* for the re¬ 
gion involved in the event (considering i?* to be the radius 
of (T Ori E), which is reasonable in the context of MOWS 
(Groote and Schmitt, 2004; Mullan, 2009) but this size 
would correspond to 10-20i?* when considering the radius 
of a K-type companion, too large for a flare associated to 
an active coronal region in the companion (Mullan, 2009). 
In addition, while flares are typical for coronal sources and 
their luminosity is compatible with that observed, flares 
of the same magnitude as those observed occurred rarely 
for low-mass stars and their typical recurrence timescale 
would have prevented several observation of such events 
over two decades (Mullan, 2009). On the other hand, 
Bouy et al. (2009) did detect the presence of a low-mass 
companion and the analysis of the source position dur¬ 
ing the Chandra-HKC flaring episode favors the low-mass 
companion as origin of the flare (Petit et al., 2012). Fur¬ 
thermore, the magnetospheric mass derived by Townsend 
et al. (2013) appears to be two orders of magnitude below 
the level required for centrifugal breakouts to happen as 
initially suggested by ud-Doula et al. (2006) . A rather sim¬ 
ilar situation was found for the magnetic Ap star IQ Aur, 
where the estimated size of the flaring region, the position 
of the X-ray source, the soft character of the recorded X- 
rays (compared to young active stars), and the high X-ray 
to radio luminosity ratio favor an origin of the flare in the 
Ap star, but the influence of a low-mass companion can¬ 
not be formally excluded (Robrade and Schmitt, 2011). 
Finally, the detection of several flares in HD 47777 was 
also reported by Naze et al. (2014) but again, it remains 
to be confirmed that they are actually associated with the 
magnetically confined winds of the star, rather than with 
an yet unknown low-mass companion of the star. 

3.3.2. Rotational modulation 

A second type of variation is rotational modulation, 
expected to occur on longer timescales than flaring. This 
type of variability strongly depends on the magnetospheric 
geometry and the orientation of the star with respect to the 
line of sight. Observations have been analyzed to reveal 
such changes. 

Amongst magnetic OB stars, rSco displays a rather 
complex field geometry, with one model restricted regions 
of densely concentrated closed loops asymmetrically dis¬ 
tributed on the stellar surface (Donati et al., 2006). Strong 
phase-locked variability could thus be expected (40% mod¬ 
ulation, Donati et al., 2006), as some regions appear in 
plain view at some phases and are occulted at others, but 
a dedicated Suzaku monitoring found at best very low am¬ 
plitude changes in flux (a few percents and mostly in soft 
band, Ignace et al., 2010). This may be explained if the 
distribution of hot gas is different (e.g. in smaller loops) 


from the magnetospheric geometry derived from optical 
spectropolarimetry. Certainly more detailed investigation 
is needed before this observational result can fully be un¬ 
derstood. 

With a simpler, mostly dipolar magnetic field, other 
magnetic OB stars should show changes related to the 
viewing angle on the confined winds (i.e. alternative edge- 
on/face-on view). Indeed, X-ray variations simultaneous 
with optical changes have been detected in at least four 
well-documented cases: 9^ Ori C (Gagne et al., 2005a; Stelzer 
et al., 2005, and references therein), HD 191612 (Naze et al., 
2010), CPD -28°2561 (Naze et ah, 2015), and NGC1624-2 
(Petit et al., 2015). For these objects, the X-ray emission 
and the Ha emission appear brighter when the confined 
winds are seen face-on. 

In the X-ray domain, variation amplitudes {Fmax /Fmin— 

1) typically reach 40-100% for 9^ Ori C , HD 191612, NGG1624- 
2, and GPD —28°2561 as well as NU Ori (Naze et al., 2014) 
and Trl6-22 (Naze et al., 2014), but the amplitude appears 
smaller in a few other cases. For example, outside flares, 
the X-ray lightcurve of crOriE shows limited variability: 
Skinner et al. (2008, see in particular their Fig. 5) de¬ 
tected variations with an amplitude (Fraax/Frain — 1) 

30% but at the 1.4% significance level only - the ampli¬ 
tude is in fact comparable to the noise and requires con¬ 
firmation. Similarly, the X-ray flux of /3 Gep changes by 
at most 10% and its phase-locked nature still need to be 
established (Favata et al., 2009; Naze et ah, 2014). Finally, 
a stable X-ray emission is observed for HD 148937, which 
is not surprising: it is always seen nearly pole-on, hence 
the view on the confined winds of this oblique rotator does 
not vary much, explaining the low amplitude changes in 
the optical domain and the absence of detectable varia¬ 
tions in X-rays (Naze et al., 2014). 

The phase-locked variability observed in X-rays could 
be explained in two different ways. On the one hand, 
it could be related to absorption. For example, if X- 
rays arise in the dense equatorial regions, then some ab¬ 
sorption of that emission could take place, reducing the 
observed (soft) X-ray flux especially when these regions 
are seen edge-on. However, to produce a large flux de¬ 
crease, this absorption increase would become noticeable 
in spectral fitting. This is the case of the extremely mag¬ 
netic 0-star NGC1624-2, where the very dense magneto¬ 
sphere produces a strong and variable additional absorp¬ 
tion (1 — 4 X 10^^cm“^, Petit et al. 2015). The spectrum 
with minimum emission clearly lacks soft X-ray photons 
when compared to that of the maximum emission phase, 
and occultation effects (see below) could only be respon¬ 
sible at most for a 20% variation in (intrinsic) flux. How¬ 
ever, this situation is exceptional, as the local, additional 
absorption due to confined winds generally remains low, as 
mentioned before, and no significant increase in absorption 
for the edge-on configuration was reported in other stars 
(Gagne et al., 2005a; Naze et al., 2014). 

On the other hand, flux variations in magnetic stars 
could be related to regular occultation of the emission re- 
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gions by the stellar body. This is the favored explanation 
and it could explain the scatter in amplitude change. For 
example, the small variations of /3 Cep are compatible with 
occultation of hot plasma located at 4-6i?*, the position 
derived from high-resolution spectra (Favata et ah, 2009). 
crOriE has a very large magnetosphere (about 31R* for 
the Alfven radius), which minimizes the impact of occul- 
tations and could contribute to explain the lack of clear 
variation detection (Naze et ah, 2014). However, as the 
flux changes are often accompanied by changes in spectral 
shape, this would imply that occultation affects plasma 
of different temperatures in different ways as could be ex¬ 
pected, e.g., if thermal stratification exists (Naze et ah, 
2014). In this context, it may be worth noting the variety 
of behaviours, which remains to be explained: HD 191612, 
Trl6-22, CPD—28°2561, and NUOri appear harder when 
brighter, while 9^ Ori C appears to be softer while brighter 
(Naze et ah, 2014). 

To assess the occultation explanation, we can model 
the confined winds by a simple optically-thin structure 
near the star, and predict the occultation degree for var¬ 
ious situations. We tested several geometries: a disk up 
to Alfven radius Ra and a ring from Ra — R* to Ra, 
which better reproduce the geometry of emitting regions 
observed in MHD simulations. Figure 5 shows how occul¬ 
tation depends on geometry (a disk is more occulted than 
a ring) and its variation with phase depends on orientation 
(i.e. for various inclination and obliquity combinations). 
Most important is the location of the plasma, of course. 
Considering a ring geometry in agreement with MHD mod¬ 
els yields a relative flux variation Emin/Fmax of ^ 0.83, or 
[Flnax/Fmin] — 1 ^ 20%, for the case of 9^ OriC while the 
observed value is [F)nax/F)nin] — 1 ^ 50%. Similar conclu¬ 
sions are reached for CPD —28°2561 or HD 191912. To get 
higher occultation degrees, one could change the plasma 
location, as structures closer to the star are more occulted 
than those further away. To get the observed value of 
6 *^ Ori C, one would however need to consider a ring of neg¬ 
ligible radius located at 1.15i?*, less than half the Alfven 
radius but also much less than the plasma location derived 
by He-like triplets (Gagne et ah, 2005b), which seems dubi¬ 
ous. For cases like Trl6-22, where [Emax/Fmin] —1 ~ 100%, 
then the hot plasma would need to be at the photosphere. 
Clearly, the simple model is underevaluating occultation, 
most probably because of its simplicity. Considering some 
limited absorption may help, but a fully optically-thick 
case would probably be inappropriate because of its phase 
shifts. 



Figure 5: Top left: Evolution with radius of the minimum value of 
the apparent flux, relative to the unocculted case, of a ring of negligi¬ 
ble Ar lying in the magnetic equatorial plane. Note that, depending 
on orientation, structures very close to the star may remain always 
partially occulted throughout the cycle, hence the y-axis value does 
not necessarily represent the observed -Fmin/Fmax, it is only the min¬ 
imum possible. Top right: Variation of relative flux with phase, for 
different brightness functions (uniform brightness, evolution as 
or The Alfven radius, inclination, and obliquity are typical 

of 0^ OriC. Bottom left: Variation of relative flux with phase, for 
different geometries. “Disk” is a structure extending from i?* to 
Ra, “ring” an annuli extending from Ra — R* to Ra‘, “Thin” cor¬ 
responds to structures lying in the magnetic equator only, “thick” 
applies to structures with a total thickness of O.SR*, “thicker” for 
thicknesses of 0.5/?*. The Alfven radius, inclination, and obliquity 
are typical of 0^ OriC. Bottom right: Variation of relative flux with 
phase, for different combinations of inclination and obliquity but a 
similar geometry (thin ring, with Ra = 2.5/?*). This agrees well 
with observations, e.g. the single minimum of 0^ OriC vs. the two 
minima of CPD —28°2561. 


^Donati et al. (2001) have envisaged the case of confined winds be¬ 
ing fully optically thick. In that case, the X-ray flux would be maxi¬ 
mum when the confined winds are seen edge-on (i.e. when the optical 
emission are minimum) as the emission from both hemispheres would 
then be seen, while only the emission of “upper half” of the confined 
winds would be observable at other phases, that of the “lower half” 
being fully absorbed. While flux variations by a factor of two have 
been detected, it must be noted that a shift by half a cycle between 
X-ray and optical maxima has not been reported, up to now. 
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3.4- Structure of confined winds, as revealed by high-resolution 
spectra 

High-resolution spectra can yield a wealth of informa¬ 
tion. With current instrumentation, line widths and shifts 
can be evaluated with precisions down to a few tens km s“^ 
in the most favorable cases (a few hundreds kms“^ more 
typically). Furthermore, the comparison of lines from H- 
like and He-like ions and of components of fir triplets 
of He-like ions constrain the temperature and location of 
the emitting plasma. However, such measurements are 
currently only possible for the brightest X-ray sources, so 
that few magnetic massive stars have been investigated 
in this respect (r Sco - Mewe et al. 2003; Cohen et al. 
2003, 9^ Ori C - Schulz et al. 2000; Gagne et al. 2005a,b, 

HD 191612 - Naze et al. 2007, HD 148937 - Naze et al. 
2008, 2012, fi Cep - Favata et al. 2009, IQAur - Robrade 
and Schmitt 2011). 

Within noise limitations, the X-ray lines of magnetic 
massive stars were found to be symmetric, and globally 
unshifted. This agrees well with MHD models. In the case 
of 9^ Ori C, global fitting however suggests small variations 
in velocity (Gagne et al., 2005a): from —75kms“^ when 
the star is seen pole-on to about 100kms“^ when seen 
edge-on. This change needs to be confirmed as the errors 
are large but also because one cannot exclude a stochas¬ 
tic variation when only a single observation per phase is 
available. If further observations provide evidence that ve¬ 
locity varies with phase, then refinement of models will be 
needed, as no such changes are currently predicted (Gagne 
et al., 2005a). 

Reported widths of X-ray lines largely depend on ob¬ 
ject and ion considered. The narrowest widths, so far, 
were found for jd Cep, whose lines are dominated by instru¬ 
mental resolution, yielding only upper limit on intrinsic 
widths (<600kms“^, Favata et al., 2009). Larger widths, 
FWHM ^ 600 — 800kms“^, were reported for ions with 
high ionization potential (Mg, Si, S) in r Sco, 9^ Ori C, and 
HD 148937, three stars with more rapid winds than /? Cep. 
Such widths are much smaller than observed for “normal” 
0-type stars {FWHM ^ Voo), where lines arise in embed¬ 
ded wind shocks distributed all over the wind, hence cover 
a larger velocity range. They indicate formation in slowly- 
moving plasma, in agreement with the confined winds sce¬ 
nario. However, most MHD models predict even narrower 
lines (Gagne et al., 2005a). 

Furthermore, lines from ions with lower ionization po¬ 
tential, notably oxygen, appear broader {FWHM ~ 1800— 
2000kms“^ Gagne et al., 2005a; Naze et al., 2007, 2008). 
These lines are associated with cooler plasma, which could 
have a different origin than hotter plasma. For example, 
the dominant hot plasma in 9^ Ori C is thought to arise 
in confined winds while the cooler one could arise in em¬ 
bedded wind shocks as in normal 0-stars (Gagne et al., 
2005a). This dual origin could be supported by the differ¬ 
ent temperatures derived from the different ions (Schulz 
et al., 2000). In Of?p stars, however, the spectra are dom¬ 
inated by the cooler component, i.e. confined winds emit 


soft X-rays in these objects (see above), but it cannot be 
excluded that current errors, which are large, are some¬ 
what blurring the picture. 

In He-like triplets, the forbidden line is suppressed when 
the density is high or the UV radiation is intense. In the 
case of massive stars, the latter effect is the most important 
one and, thanks to dilution with distance, enables us to 
locate the emitting region. In t Sco, 9^ Ori G, HD 148937, 
and fd Cep, the start of the emitting region is found to 
be close to the photosphere, at radii r ~ 1.5 — 3 i?* for 
the first three stars and r ^ 4 — 6 R* for the latter case. 
These values are only slightly lower than the corresponding 
Alfven radii of these stars, they thus appear qualitatively 
compatible with MHD simulations. In IQ Aur, the forbid¬ 
den line was found to be normal, suggesting a formation 
radius larger than 7i?* (Robrade and Schmitt, 2011) - de¬ 
spite a large value, this also appears compatible with the 
supposed location of confined winds in this star. 

4. Discussion, Conclusion and Future Outlook 

The discovery of magnetism in massive O stars is a 
relatively recent phenomenon although magnetism in B 
stars was known for much longer. This led to a number of 
new observational campaigns such as MiMeS (Wade et al., 
2012) or BOB (Morel et al., 2015) that increased the total 
number of known magnetic OB stars significantly. There 
have also been some theoretical developments in an at¬ 
tempt to understand this increasingly larger group of stars. 
In this brief review, we discussed specifically the X-ray 
emission from these objects. 

Although there are some outliers, semi-analytic scaling 
of XADM paradigm can correctly estimate X-ray luminosi¬ 
ties from most magnetic OB stars. A phenomenon called 
shock retreat can reduce the X-ray luminosity, and for the 
case of B stars with extremely low mass-loss rate, it may 
even quench the wind within the closed magnetosphere. 
Rotation can counter some of the effects, but much work 
remains to be done in order to fully understand these com¬ 
plex objects. 

As the number of magnetic detections increases with 
time, future X-ray missions will certainly enlarge the sam¬ 
ples. Indeed, most of the results presented here come from 
a relatively small size sample while a large number of pa¬ 
rameters are involved in the confined wind phenomenon 
and its observations (magnetic strength, magnetic geome¬ 
try, inclination of the system, and stellar parameters such 
as temperature, rotation rate, wind velocity and mass-loss 
rate). Notably, only about 10 magnetic 0-stars are cur¬ 
rently known and have been studied in X-rays. Moreover, 
the sensitivity of current facilities does not allow a large 
study of magnetic massive stars at high spectral resolu¬ 
tion, nor the detection of very faint X-ray emission (e.g. 
Ap/Bp stars). This certainly poses a challenge to under¬ 
stand these object fully and with confidence. 

The advent of the new X-ray mission Athena-f (Nan- 
dra et al., 2013) may notably enable us to study con- 
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fined winds in radically different environments like low- 
metallicity galaxies such as the Magellanic Clouds, where 
massive stars are known to display weaker stellar winds. 
Such new data would allow to test our current models 
on vastly different cases, thus improving their reliability. 
The very first extragalactic magnetic massive candidates 
have already been proposed (Maitzen et ah, 2001; Paun- 
zen et ah, 2011; Naze et ah, 2015; Walborn et ah, 2015), 
so that several targets will become available in the next 
decade for such studies. 

Another important observational information comes from 
monitoring the X-ray emission. It is a natural way to probe 
the magnetosphere, and it provides strong constraints on 
its hot plasma content. Unfortunately, such programs are 
difficult to obtain because of their length and the need of 
time-constraints, therefore only few objects could be fol¬ 
lowed in detail. Moreover, when a monitoring exists, it 
generally covers a single cycle or less - it is thus currently 
difficult to separate potential stochastic changes from pe¬ 
riodic ones. With a better sensitivity, such as foreseen for 
Athena-|-, shorter exposures will be needed, reducing the 
total duration hence facilitating the gathering of multiple 
exposures. 

Finally, high spectral resolution data also provide cru¬ 
cial information by allowing detailed studies of the X-rays 
associated to confined winds. Currently, only a few such 
spectra were obtained, and they did not yet allow analy¬ 
ses of the line profiles, as lines are barely resolved. Future 
facilities such as the European Space Agency’s Athena-I- 
observatory (Nandra et ah, 2013) should provide access 
not only to high spectral resolution for more objects but 
also to better spectral resolution than currently available. 
This will allow us to pinpoint the properties of a larger 
sample with better quality, further constraining the hot 
plasma distribution and kinematics. Of utmost interest 
in this context are monitoring performed at high resolu¬ 
tion. The behaviour of different lines will then be com¬ 
pared, e.g. those from ions with low-ionization potential 
and those associated with high ionization potential. Any 
difference between them will indicate different emission 
processes, which will be identified thanks to detailed line 
profile studies. Furthermore, if velocity shifts such as ten¬ 
tatively identified in 6^ OriC are confirmed, then detailed 
tomographic analyses will become possible, enabling the 
reconstruction of the magnetospheric activity kinematics. 

Along with the observational developments, theory must 
also keep up the pace. Currently, most of the fully consis¬ 
tent MHD simulations are limited to 2D. One exception is 
the 3D MHD study of Ha emission from 9^ Ori C by ud- 
Doula et al. (2013). Since rotation in this case was dynam¬ 
ically unimportant, rotation and field axes were assumed 
to be aligned. However, in many other cases rotation is 
dynamically important and the field is tilted with respect 
to the rotation axis. Such oblique rotator cases must be 
modelled using full 3D MHD, something currently missing. 
Moreover, for stars for which magnetic confinement pa¬ 
rameter 77* ^ 1, MHD cannot be used due to high Alfven 


speed that forces Curant time to be impractically small, 
rendering such numerical simulations essentially impossi¬ 
ble. 

In such extreme cases, alternate methods must be em¬ 
ployed. Townsend and Owocki (2005) present a semi- 
analytic method for modelling the circumstellar environ¬ 
ment of early-type stars, which is essentially a generalized 
extension of Babel and Montmerle model. They assume 
that the magnetic field lines remain rigid and corotate with 
the star without being influenced by the dynamics of the 
wind. This is essentially equivalent to the assumption of 
77* —>■ 00. This Rigidly Rotating Magnetosphere (RRM) 
model can be applied to an arbitrary field geometry and 
tilt angle between the field and rotation axes. However, 
there are no dynamic forces involved in this model, instead 
effective gravitational -I- centrifugal potential is calculated 
based on the constrained motion of the plasma. Accu¬ 
mulation surfaces are assumed to lie along the location of 
minimum effective potential. 

As a further improvement to the RRM model, Townsend 
et al. (2007) introduced a new Rigid-Field Hydrodynamics 
(RFHD) method to modelling the circumstellar environ¬ 
ments of strongly magnetic massive stars as defined above. 
Just like in the RRM method, the field lines are treated 
as rigid, but now the flow along the lines is computed 
self-consistently using hydrodynamical equations includ¬ 
ing line force due to radiation. In this ansatz, flow along 
each each field line is considered to be an independent ID 
flow. They perform a large number of such ID calculations 
for differing field lines, then piece them together to build 
up a time-dependent 3D model of a magnetosphere. 

This method is now superseded by Arbitrary Rigid- 
Field Hydrodynamics (ARFHD) method which now allows 
any configuration of self-consistent magnetic topology. It 
also improves some of the numerical algorithms and radia¬ 
tive cooling and includes thermal conduction (Bard and 
Townsend, 2015). Since the flow along each field line can 
be solved independently of other field lines, the computa¬ 
tional cost of this approach is a fraction of an equivalent 
full MHD simulation. However, this method suppresses 
any coupling between the flow on adjacent field lines which 
in practice may have some important ramifications. As 
such, a new, improved method will be needed in the fu¬ 
ture. 
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